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Abstract 

The phenomenon of spontaneous symmetry breaking admits a physical interpre- 
tation in terms of the Bose-condensation process of elementary spinless quanta. In 
a cutoff theory, this leads to a picture of the vacuum as a condensed medium whose 
excitations might deviate from exact Lorentz covariance in both the ultraviolet and 
infrared regions. For this reason, the conventional singlet Higgs boson, the shifted 
field of spontaneous symmetry breaking, rather than being a purely massive field, 
might possess a gapless branch describing the long-wavelength fluctuations of the 
scalar condensate. To test this idea, that might have substantial phenomenological 
implications, I compare with a detailed lattice simulation of the broken phase in the 
4D Ising limit of the theory. The results are the following: i) differently from the 
symmetric phase, the single-particle energy spectrum is not reproduced by the stan- 
dard massive form ii) for the value of the hopping parameter k = 0.076, increasing 
the lattice size from 20^ to 32^, the mass gap is found to decrease from the value 
0.392(1) reported by Balog et al. (see Nucl. Phys. B714 (2005) 256) to the value 
0.366(5). Both results confirm that, in the infrared region, the standard singlet 
Higgs cannot be considered as a simple massive field. Several arguments indicate 
that, approaching the continuum limit of the lattice theory, the observed volume 
dependence of the mass gap might require larger and larger lattice sizes before to 
show up. 



1. Introduction 



The idea of a 'condensed vacuum' is generally accepted in modern elementary particle 
physics. Indeed, in many different contexts, one introduces a set of elementary quanta 
whose perturbative empty vacuum state |o) is not the true ground state of the theory. In 
the physically relevant case of the Standard Model, where the condensation process cor- 
responds to spontaneous symmetry breaking (SSB) through the A$'^ sector, the situation 
can be summarized saying [1] that "What we experience as empty space is nothing but 
the configuration of the Higgs field that has the lowest possible energy. If we move from 
field jargon to particle jargon, this means that empty space is actually filled with Higgs 
particles. They have Bose condensed" . The translation from field jargon to particle jargon 
can be obtained, for instance, along the lines of Ref.[2] where the substantial equivalence 
between the one-loop effective potential of quantum field theory and the energy density 
of a dilute particle gas was established. 

With these premises, it was argued in Ref. [3] that such a physical representation 
of the vacuum might run into a contradiction with the conventional picture where the 
standard singlet Higgs boson, the shifted field of SSB, is represented as a purely massive 
field. In fact, in the hydrodynamical regime of disturbances whose wavelengths are much 
larger than the mean free path of the elementary constituents, quite independently of 
the Goldstone phenomenon and even for a spontaneously broken one-component 
theory, one would expect that the lowest excitations of the system correspond to small 
displacements of the condensed quanta that already preexist in the ground state. In this 
case, for momenta p — > 0, the energy spectrum should end with an acoustic "phonon" 
branch ~ Cs|p| as it happens in superfiuids and in all known condensed media. The 
existence of such a regime is, in fact, a very general result [4] that does not depend 
on the details of the short-range interaction and even on the nature of the elementary 
constituents. For instance, the same coarse-grained description is found in supcrfiuid 
fermionic vacua [5] that, as compared to the Higgs vacuum, bear the same relation as 
superfiuid ^He has to superfluid ^He. 

In this sense, as discussed in Ref. [6], the idea of a massive excitation spectrum of the 
form a/p^ + down to p = can be considered equivalent to the incompressibility limit 
where there is no phonon branch and = +oo. This can be a valid representation of the 
vacuum for the continuum theory where, as in axiomatic quantum field theory, one can 
derive the exact Lorentz covariance of the energy spectrum. However, in the condensed 
phase of a cutoff theory one is faced with "reentrant violations of special relativity in the 
low-energy corner" [7] that only disappear in the continuum limit where the ultraviolet 
cutoff A ^ oo. 

In a condensate of spinless particles, these violations would be simple density fiuc- 
tuations extending over a small region |p| < 5. Lorentz covariance can become exact 
if the shell 5 vanishes, in units of the scale Mh associated with the massive part of the 
spectrum, in the limit A — oo. Introducing dimcnsionless units, e = -g^ and t = the 
"reentrant" nature of the Lorentz violating infrared effects implies the equivalence of the 
e — > and t — > oo limits as, for instance, embodied in the single relation [6] 



In the cutoff theory, the very long- wavelength phonons of the Bose condensate would 
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produce a very weak 1/r potential [8]. Thus one can consider [6] possible viable phe- 
nomenological frameworks to constrain the relative magnitude of A and S, while keeping 
M/i at the Fermi scale. For instance, for = 250 GeV and A ~ 10^^ GeV one gets 
S ~ 10~^ eV and a ratio e = S/Mh ~ 4 • 10'"^''' that might well represent the physical 
realization of a formally infinitesimal quantity. If this were the right order of magnitude, 
the non-Lorentz covariant density fluctuations of the vacuum would start to show up 
from wavelengths larger than a few millimeters up to infinity. These lengths are actually 
infinitely larger than the Fermi scale but, nevertheless, have a physical meaning. 

Now, one might ask if there is any "first- principle" support for the intuitive physical 
scenario sketched above. To this end, two arguments were given in Ref.[3] to motivate 
a possible failure of the standard perturbative picture in the far infrared region of the 
broken-symmetry phase. The first argument starts from the separation of the full scalar 
field into a constant background (p and a shifted fiuctuation field h{x) 

<^{x)^(f> + h{x) (2) 

Here, in order the separation to be unambiguous, (p denotes the average field in a large 
four-volume Q, i.e. 

d'x (3) 

and the limit — > oo has to be taken at the end. In this way, the full functional measure 
can be expressed as 

r r+oo r 

/ [ri$(x)] = / d(t) [dh{x)] (4) 



where the functional integration on the right-hand side is over all quantum modes with 
four-momentum ^ 0. 

In the standard approach, one describes the shifted field h{x) as a purely massive 
field whose mass is related to the quadratic shape of a non-convex semiclassical effective 
potential Vnc(0) at one of its absolute minima. This identification is based on introducing, 
for any given 0, the zero- momentum two-point function T2{p = 0), i.e. the inverse zero- 
momentum connected propagator D^^{p = 0), through the relation 

Z,-(p=0,^^ (5) 

and defining M| from the value of D^^{p = 0) at = ±f , the absolute minima of Hvfc(0)- 
In this series of steps one treats 0, the zero-momentum mode of the full scalar field 
as a purely classical object without taking into account that, on the base of Eq.(4), 
there should be one more functional integration over d(f). In an alternative approach, one 
can start from the generating functional in the presence of a constant source J 

/ + 00 
d(t> exp[-l](yNc(0) - J4>)] (6) 
-oo 

where Vnc(0) is the non-convex effective potential obtained, order by order in the loop 
expansion, after integrating out all modes with 7^ 0. 

Further, introducing the generating functional for connected Green's functions w{J) 
through 

li»W=ln|^ (7) 
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the vacuum field and the zero-momentum propagator are defined from the relations 

dw , , 

'AS) = 5J (8) 

and 

G(p=0) = ;j^ (9) 

in the double limit J — > and f2 — > oo. 

Now, it is easy to show that in order to have SSB, i.e. 

^lim^ (^(J) 7^ (10) 

the four- volume Q has to diverge in such a way that the dimensionless quantity x — Q.Jv 
tends to a non-zero limit. However, even assuming that x tends to ±cxd, so that the vacuum 
field (/?(J) tend to ih?;, still the zero-momentum inverse propagator G~'^(j) = 0) does not 
uniquely tend to M^. In fact G~^{p = 0) is the second derivative of the exact effective 
potential that, as defined from the Legendre transform of the generating functional, is 
not an infinitely differentiable function in the infinite- volume limit [9]. For this reason, in 
the broken-symmetry phase, there are two solutions 

G,\p = 0) = M', G^\p = 0) = (11) 

as if there were two different particles in the theory and not just one. 

In this sense, the situation of the broken phase is reminiscent of what happens in 
superfluid ''He. Also there the spectrum is considered to arise from the combined effect 
of two different types of excitations, phonons and rotons, whose separate energy spectra 
match giving rise to a complicated pattern. Just following this analogy, one can try to 
use quantum hydrodynamics [10] to obtain a physical interpretation of the parameter 
associated with the massive part of the spectrum in terms of the energy-gap for vortex 
formation in a superfluid medium possessing the same constituents and the same density 
as the scalar condensate. In this picture, where phonons would be physically cut off above 
those infinitesimal momenta that correspond to the inverse-millimeter range, the branch 
of the spectrum ~ ^/p^ + covers in practice the full range that is relevant for particle 
physics. 

The above derivation of two solutions for the zero-momentum propagator is non- 
perturbative and independent of any diagrammatic expansion. In Ref.[3], however, it was 
suggested that the same conclusion is also obtained taking into account the effect of the 
one-particle reducible, zero-momentum tadpole graphs. These enter the diagrammatic 
expansion for the propagator in the presence of a constant background field and can be 
considered a manifestation of the quantum nature of the scalar condensate. These effects 
are usually neglected in the standard procedure where, at the minima of Vnc(0), one 
defines the connected propagator from a Dyson sum of IPI graphs alone. The reason is 
that the one-particle reducible tadpole graphs are proportional to the one-point function 

J(^) ^ ^ (12) 

that vanishes at (f) — ±v. However, these zero-momentum graphs are attached to the 
other parts of the diagrams through zero-momentum propagators. Therefore, ignoring 
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their contribution at (f) — ±v (or including their effect in a pure perturbative way with a 
form of the zero-momentum propagator G(p = 0) = + ...) is equivalent to assume the 

h 

regularity of G{p = 0) when J ^ 0. Instead, if one relaxes this assumption and looks for 
the general form of G~^{p = 0), one finds again two solutions as in Eq.(ll). 

Addressing to Ref.[3] for more details, the tacit assumption at the base of the standard 
approach is better illustrated with a very simple example. Consider the quadratic equation 

f-\x) = l + x^~g^x^f{x) (13) 

for (yf^ ^ 1. The analogy with the problem of the one-particle reducible zero-momentum 
tadpole graphs is established when comparing f{x) with G{p = 0), at a given value of 0, 
and the limit x — > with the limit J(0) 0. Standard perturbation theory is based on 
the iterative structure /reg = 1/(1 + x^) + O^g"^) that provides a class of solutions that 
are regular for a; — where /rcg(O) = 1. For this class of solutions, that corresponds to 
the massive propagator as defined from the one-particle irreducible graphs only, the third 
term in the r.h.s. of Eq.(13) vanishes identically for x — > 0. On the other hand, Eq.(13) 
has also a singular solution fgi^g ~ 1/ g^x^ for x — > and this corresponds to a divergent 
zero-momentum propagator when (j) ±v. This can only be discovered by retaining the 
full non-linearity of the problem where the zero-momentum propagators joining to the 
vacuum sources are not approximated perturbatively 

I realize that the various arguments given above cannot be considered a "proof that 
there is a gapless branch in the excitation spectrum of the broken phase. They represent, 
however, convergent indications that the standard perturbative derivation of a massive 
spectrum, which is based on the quadratic part of the shifted classical lagrangian after 
simply replacing = ±v in Eq.(2), may be too naive. A deeper understanding of the 
ground state might be needed to obtain a proper description of the far infrared region 
of the condensed phase. For instance, exploiting the superfiuid analogy, the methods of 
quantum hydrodynamics (sec e.g. Ref.[ll]) might represent a natural alternative. 

After this general introduction, 1 will provide in the rest of the paper two more argu- 
ments that also indicate an unconventional infrared behaviour. The first new argument, 
presented in Sect. 2, is based on the results of Ref.[12]. It uses the numerical solution of 
the coupled RG equations for the effective potential and field strength at various values 
of the infrared cutoff k. Expressing the full scalar field as in Eq.(2), the numerical results 
indicate that, in the broken phase, the fluctuation field h{x) cannot be represented in 
terms of purely massive states in the k limit. 

The second new argument is based on a numerical simulation of the broken phase 
in the Ising limit. If the shifted fluctuation fleld were a purely massive fleld, the single- 
particle energy spectrum should be reproduced by (the lattice version of) the standard 
massive form -|- const, in the p — limit. At the same time, there should be no 
observable change in the mass-gap increasing the linear lattice size above a typical length 
scale associated with 7-8 correlation lengths. As discussed in Sect. 3 (details in the Tables 
at the end of the paper) , both expectations are not consistent with the lattice data. 

Finally, Sect. 4 will contain a summary and the conclusions. 

2. The RG equations for \/(0) and ^(0). 

The aim of Ref.[12] was to study the effective potential V{(f)) and the fleld strength Z(0), 
as functions of the background fleld 0, at various values of the infrared cutoff. This is a 
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widely accepted technique where one starts from a bare action defined at some ultraviolet 
cutoff A and effectively integrates out shells of quantum modes down to an infrared cutoff 
k. This procedure generates a A;— dependent effective action rfc[$] that evolves into the 
full effective action r[$] in the /c — > hmit [13-17]. The A;— dependence of Tk[^] is 
determined by a differential functional flow equation that is known in the literature in 
slightly different forms. In particular, with the flows discussed in detail in Rcf. [18] one 
starts form flrst principles and obtains a class of functionals that interpolates between the 
classical bare Euclidean action and the full effective action of the theory. However, some 
features, such as the basic convexity property of the effective action for A; — > [19-22], 
are independent of the particular scheme. 

In this approach, the relevant quantities are the A;— dependent effective potential Vk{4>) 
and fleld strength Zk{(l)), which naturally appear in a derivative expansion of Tk[^] around 
a space-time constant conflguration $(a;) = 4>. They are governed by two coupled equa- 
tions derived and discussed in Refs. [22-25]. 

Introducing dimensionless variables : t — ln(A/A;), x — k^'^^'^cf), V{t,x) — k~^Vk{(t>) 
and Z{t,x) = Zk{(J)) (where D indicates the number of space-time dimensions) and defln- 
ing the flrst derivative of the effective potential f{x,t) = dxV{t,x), these coupled equa- 
tions can be expressed in the form [12] (/' = d.,:f{t,x), Z' = dxZ{t,x),...) 

dl ^ {D + 2) {2-D) dl _ ^_ ^ 

dt 2 2 dx (47r)^/2 dx ^ ' 

dZ_ _ {2-D) ^dZ^ 1 dfZ'^_,„ 



dt 2 ' dx (47r)^/2 dx\Z 

-f/z ( f{^z'f 18D - ^2 - 20 {A-D)Z'f" (/ 



(47r)^/2 \^ Z3 ^ 24 Z^^ 6Z2 ~QZ^ j ^^^^ 

It is easy to show that these two coupled equations can be transformed into the structure 
(^,J=l-3) 

where the components of the vector Ui{x,t) are the unknown functions of the problem 
and where Pjj, Qi and Ri can depend on x, t, Ui, In this way, the numerical solution 
was obtained with the help of the NAG routines. 

The analysis of Ref.[12] was focused on the quantum-fleld theoretical case D = 4 
assuming standard boundary conditions at the cutoff scale: i) a renormalizable form for 
the bare, broken-phase potential 

Vj,{<P) ^ ~M'<P' + X<f>' (17) 

and ii) a unit normalization condition for the derivative term in the bare action 

Za(0) = 1 (18) 

It was also assumed a weak-coupling limit A = 0.1, flxing M = 1 and A = 10. In this 
way, there is a well deflned hierarchy of scales where the infrared region corresponds to 
the limit A; < M < A. 
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The numerical results can be summarized as follows. For not too small values of the 
infrared cutoff k, the effective potential Vk{4>) remains a smooth, non-convex function of 
4> as in the loop expansion. In this region of k one also finds a field strength Zk{(j)) ~ 1 
for all values of (p. 

However, a tiny scale S ~ 0.15 exists such that ior k < S the effective potential Vk{(f)) 
starts to flatten in an inner region of \(f)\ while still matching with an outer, asymptotic 
shape of the type expected in perturbation theory. The flattening in the inner |0|-region, 
while reproducing the expected convexity property of the effective potential, does not 
correspond to a smooth behaviour. 

For such small values of k there are large departures of Zk{(f)) from unity in the inner 
|0|— region with a strong peaking at the end point |0| = |0(/c)| of the flattening region. 
On the base of the general convexification property, the A; — > limit of such end point, 
0(0), coincides with one of the minima ±v of a suitable semiclassical, non-convex effective 
potential and is usually taken as the physical realization of the broken phase. 

To interpret these results, let us start from the usual point of view where SSB is 
described in terms of a classical potential with perturbative quantum corrections. These 
corrections, with the choice of the bare parameters of Ref.[12], are typically small for 
all quantities. In particular Z, in perturbation theory, is a non-leading quantity since 
its one-loop correction is ultraviolet finite. Therefore, perturbation theory predicts tiny 
deviations of Z from unity, ~ 10~^. This is also in agreement with the assumed exact 
"triviality" of the theory [26] that requires Z ^ 1 m the continuum limit. 

This prediction fits well with the profile of Zk{(/)) for not too small values of the infrared 
cutoff. However, for k < 6, there are large deviations from unity with the mentioned strong 
peaking phenomenon. 

If we express the full scalar field as in Eq.(2), the above results indicate that 

the higher frequency components of the fluctuation fleld h{x), those with 4- momentum 
Pfj_ such that S < \p\ < A, represent genuine quantum corrections for all values of the 
background field 4> in agreement with their perturbative representation as weakly coupled 
massive states. 

On the other hand, the components with a 4-momentum such that \p\ < 5, are non- 
perturbative for values of the background field in the range < (f) < (f){\p\) . In particular, 

the very low-frequency modes with \p\ behave non-perturbatively for all values of 
the background in the full range < |0| < f and thus cannot be represented as standard 
massive states. In fact, non-perturbative infrared phenomena cannot occur in a genuine 
massive theory. 

Notice that the unexpected effects show up in connection with the convexification 
process, precisely as discussed in Ref.[3] and reviewed in the Introduction. In this sense, 
one can say that the convexification process is induced by the infinitely long-wavelength 
modes that, so to speak, "live" in the full region < |0| < v. 

Notice also that, by itself, the existence of a non-perturbative infrared sector in a 
region < \p\ < S might not be in contradiction with the assumed exact "triviality" 
property of the theory if, in the continuum limit, the infrared scale S vanishes in units 
of the physical parameter Mh associated with the massive part of the spectrum. This 
means to establish a hierarchy of scales 6 -C <^ A such that ^ when ^ ^ 0, 
as discussed in the Introduction. 

If this happens, the region < \p\ < 5 would just shrink to the zero-measure set 
Pf^ — 0, for the continuum theory where Mh sets the unit mass scale, thus recovering the 
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exact Lorentz covariance of the energy spectrum since the point = forms a Lorentz- 
invariant subset. In this hmit, the RG function Zk{4>) would become a step function 
which is unity for all finite values of k (and 0) and is only singular for A; = in the range 
< 101 <v. In this way, one is left with a massive, free-field theory for all non-zero values 
of the momentum, and the only remnant of the non-trivial infrared sector is the singular 
re-scahng of (the projection of the full scalar field $(x) onto — 

3. Lattice simulation of the broken phase in the Ising limit 

In this section I will present the results of a numerical simulation of the theory. The main 
goal of the numerical experiment is to check the standard assumption that for k — > 
the excitation spectrum of the broken phase is the same as in a weakly coupled massive 
theory. Quite independently of the theoretical arguments given in the Introduction and 
in Sect. 2, the idea of deviations from the simple massive behaviour of perturbation theory 
is supported by the lattice results of Ref.[27]. There, differently from what happens in 
the symmetric phase, the connected scalar propagator was found to deviate significantly 
from (the lattice version of) the massive single-particle form l/(p^ + const) for 0. In 

particular, looking at Figs. 7, 8 and 9 of Ref.[28], one can clearly see that, approaching the 
continuum limit of the lattice theory, these deviations become more and more pronounced 
but also confined to a smaller and smaller region of momenta near p^ — 0. 

As in Ref.[27], the test has been performed in the Ising limit that traditionally has 
been chosen as a convenient laboratory for the numerical analysis of the theory. In this 
limit, a one-component $| theory becomes governed by the lattice action 

'S'ising = XI XI [^(^ + ^n)(P{^) + - ^n)(P{^)] (19) 

X [X 

where 0(x) takes only the values ±1. The broken-symmetry phase corresponds to values 
of the hopping parameter k> Kc with ~ 0.07484 [29]. 
Observables include the bare magnetization: 

VB^m) , 4>^Y.^{x)lL'' (20) 

X 

(where is the average field for each lattice configuration) and the bare zero-momentum 
susceptibility: 

xiatt = LM<H'>-(H)T- (21) 

The equivalence of these two definitions with other standard lattice definitions that are 
found in the literature has been discussed in Ref . [30] . 

Let us now consider the mass gap of the theory, say ttits (k = 0) , that is traditionally 
extracted at zero 3-momentum from the exponential decay (TS='Time Slice') of the 
connected two-point correlator. This is a well known strategy (see for instance Ref. [31]) 
that, however, for the sake of clarity it might be useful to review in some detail. 

The Fourier transform of the connected two-point correlator can be expressed as 

Ci(t,0;k) = (S'e(t;k)5e(0;k) + 5,(t;k)5,(0;k))eonn (22) 

where 

Scit: k) = ^ X '^(^' ^) co«(^^ • ^)' (23) 
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5,(i;k) = -l J]</.(x,i)sin(k.x). (24) 

X 

Here, < t < is the Euclidean time; x is the spatial part of the site 4- vector x^; k is 
the lattice 3-momentum k = {2-^/ L){nx,ny,nz), with {nx,ny,nz) non-negative integers; 
and (...)conn denotes the connected expectation value with respect to the lattice action, 
Eq. (19). 

To obtain the time-slice mass, one starts from the general expression (see Ref. [31]) 

Ci(t,0;k) = 5^|^,|V^««* (25) 

a 

where the sum is over the eigenstates of the lattice Hamiltonian corresponding to the given 
value of k. In perturbation theory, these are represented in the Fock space as massive 
a— particle states coupled to total momentum k and one predicts -E'i(k) < -E'2(k) < .... 
Therefore, using the lattice dispersion relation 



(k) = 2(cosh£;i(k) - 1) - 2^ (1-cosA;^) (26) 



one can extract the mass from the single-particle energy -E'i(k). 

Now, in an interacting theory, a model-independent approach to the energy spectrum 
would require to extract the leading term Ei (k) from the exponential decay at asymptotic 
t. However, in a real numerical simulation the vanishing of the correlator can hardly be 
studied for asymptotic times with enough statistics. Thus, in order to extract the leading 
single-particle energy, one is forced to introduce some model-dependent assumptions, such 
as the shape of £'i(k), its relation with the higher excitations -E2(k), £'3(k), .. and so on. 

On the other hand, in a weakly coupled theory (as in a "trivial" theory close to the 
continuum limit), there is a simple strategy to extract the single particle energy that has 
the advantage of being free of uncontrolled theoretical assumptions. It can be applied 
once the ratio 

is expected to be very close to unity. In this regime, by further noticing that the residual 
corrections to a single-particle correlator are of order 

(1 - i?)e-^"('^)* (28) 

with A„(k) = £^n(k) — Ei{k) > and n=2,3,. . . , and thus are suppressed by additional 
phase-space factors, one can start parameterizing the full correlator (with periodic bound- 
ary conditions) in terms of an effective single-particle spectrum 

Ci{t, 0; k) = A[exp{-E{k)t) + exp(-E(k)(Lt - t)) ] (29) 

Whenever higher excited states were really needed to describe the lattice data, one should 
detect appreciable differences, both in the value of the normalized chi-square and in the 
fitted value of -E(k), by simply varying the range of t where the fit is performed. This 
means to look for the stability of the results finding a characteristic 'plateau' E^P\'k) 
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(p='plateau') for the fitted energy where the different indications converge and that can 
be used to safely extract the time-shce mass. 

After these prehminaries, let us now consider the numerical results. The simulations 
were performed with a numerical code written by Paolo Cea and Leonardo Cosmai. It 
employs the Swendsen-Wang cluster algorithm [32], using the cluster improved estima- 
tor [33] to compute the time slice correlations. The statistical analysis is performed using 
the jackknife method [34]. 

To illustrate with an example the determination of the energy plateau, let us first 
look at the symmetric phase for k — 0.074. The raw data for the connected correlator 
are reported in Tables 1, 2 and 3 for three values of the lattice 3-momentum = 0, 
k2 = 0.375 and = 0.922 and two lattice sizes. I also show the quality of the fits and 
the resulting time-slice mass in various ranges of t. As one can see, after simply skipping 
the first one or two time slices, one gets very good stability with mass values that are 
independent of the spatial momentum and in excellent agreement with each other. This 
confirms to a very high precision that, in the symmetric phase, the single-particle energy 
spectrum of the theory is well reproduced by a standard massive spectrum. 

Let us now consider the broken phase. This was simulated choosing the value of the 
hopping parameter n = 0.076. At this value of k, in fact, published data by Balog et al. 
[30] from a 20^ lattice provide the estimate mTs(O) — 0.392(1). A check that this is a real 
mass gap can be obtained in two ways. First, as done for the symmetric phase, one can 
study the stability of mTs(k) at various values of the spatial momentum. Second, one can 
check the dependence on the lattice size. In fact, if mTs(0)=0. 392(1) were a real mass 
gap, a 20^ lattice contains already ~ 8 correlation lengths and, thus, there should be no 
significant change by further increasing the lattice size. 

The simulations were first performed on 20^ lattices at different values of k^. In Table 
4, I report the results of a simulation of 6 Msweeps for k = 0. As one can check, the 
magnetization and the susceptibility are in excellent agreement with those reported by 
Balog et al. for the same lattice size (see the corresponding entries for k. = 0.076 in Table 
3 of Ref.[30]). 

Concerning the time-slice mass, the value mTs(O) =0.3920(24) is also in excellent 
agreement with the result of Ref.[30]. In particular, using all data in the t-range 1-19 the 
fit with Eq. (29) goes through all central values with a total chi-square which is less than 
0.1. This shows that, for k = 0, the deviations of the correlator data from a pure single 
exponential are confined to the first time slice. 

The results for k 7^ are reported in Tables 5-9. The number of sweeps is such to 
provide statistical errors of the time-slice mass that are approximately equal to those 
reported in Table 4 for k = 0. As one can see, differently from the symmetric phase, there 
is a distinct dependence of mTs(k) on the spatial momentum (that was already observed 
in Ref.[27]). 

Let us first compare with the data in Tables 5 and 6. These refer to the lowest 
two momenta of a 20^ lattice. Again, as for k = 0, the deviations from a pure single 
exponential are limited to the first time slice since the fit in the full t-range 1-19 with Eq. 
(29) goes through all central values to a very high precision. Using the results of the fit 
in this range one obtains time-slice masses mTs(k) = 0.3992(21) and 0.4039(21) for the 
integer assignments of the 3-momentum (1,0,0) and (1,1,0) respectively. These values are 
not consistent within 3-5 a with the corresponding value 0.3920(24) reported in Table 4. 
Therefore, the lattice data of the broken phase are not well reproduced by the standard 
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massive form in the k — > limit. 

Let us now consider the data at the higher spatial momenta reported in Tables 7-9. 
From the substantial equivalence of the fits in the ranges 2-18 and 3-17, one can deduce 
that, in this momentum range, the deviations from a pure single exponential affect now 
the first two time slices. However, notice the difference with the symmetric phase. Prom 
Table 3, using the fit results in the t-range 2-18, one obtains mTs(k) =0.2133(33) in very 
good agreement with the mass values reported in Table 1. From Tables 7-9, on the other 
hand, one obtains the corresponding entries mTs(k) = 0.4122(15), 0.4127(18), 0.4129(20) 
that are not consistent within 7-8 a with the mass value 0.3920(24) of Table 4. Using the 
fit results in the t-range 3-17, statistical errors become larger but the discrepancy remains 
at the level of 3-5 a . 

On the basis of the previous results for k ^ 0, it should be clear that interpreting 
the value mTs(O) = 0.392(1) as a real mass gap is not so obvious. Actually, if one looks 
for stability with respect to changes of the spatial momentum, it is only at higher mo- 
menta that one gets consistency with the standard concept of " mass" as a k-independent 
quantity. In fact, looking at Tables 7-9, if one requires both a good-quality fit and a sub- 
stantial independence of the chi-square per degree of freedom on the sample of data (as it 
happens for the fits 2-18 and 3-17) one obtains a remarkable momentum independence for 
a mass value ~ 0.412(3). This confirms the basic idea of the Introduction, namely that 
the shifted fiuctuation field is not purely massive and that the end point of the spectrum 
at k = cannot be used to extract the "mass". 

The idea that mTs(O) might not be a real, physical mass gap is also confirmed by its 
volume dependence. To this end, I report, in Tables 10-13 the raw correlator data from 32^ 
lattices. The total statistics is 6 Msweeps as for the simulation at k = on the 20^ lattice 
reported in Table 4. To shorten the total time of the numerical experiment, however, the 
full statistics, for the 32^ case, was collected using different computers. In the 32^ case, I 
only report data for which the S/N ratio is larger than unity. 

Notice that the magnetization and the susceptibility are completely consistent with 
those measured on the 20^ lattice. On the other hand, looking at the results of the global 
fit to the correlator data, one finds a serious discrepancy in the value of the mass gap. In 
fact both from Table 4 and Table 14, there are clean indications for the existence of an 
energy plateau. However the two values are substantially different. In particular, the 32"^ 
value £^*^^)(0) ~ 0.364(5) gives a time slice mass mTs(O) = 0.366(5) that is not consistent 
within 5 a with the value 0.392(1) reported by Balog et al. for the 20^ lattice. 

As an additional check, a fit with 2 exponentials was also performed along the lines 
indicated in Ref.[35]. In this case, Eq.(29) is replaced by a constrained 2-mass formula 
with -^2(0) = 2i?i(0). The results of this other type of fit are shown in Tables 15 and 
16. The entries that are unchanged correspond to fits where the normalization of the 
2-particle term is set to zero by the fit routine. For other entries, the statistical errors 
become much larger. However, the discrepancy between the mass value from the 20^ 
lattice and that obtained from the 32^ lattice remains. 

Notice that, differently from mTs(O) the susceptibility xiatt remains almost unchanged 
when increasing the lattice size. Since the susceptibility is nothing but the zero-four- 
momentum connected propagator, one can express the numerical stability of the result 
in the form xiatt — Ga{p — 0) thus relating xiatt to the a- type, massive solution for the 
zero-momentum propagator in Eq.(ll). On the other hand, the lattice results suggest also 
that, if there were a gapless branch of the shifted field, the lattice definition of the b-type 
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of solution might be identified through the relation m^g(O) = G'^^(j) — 0). To prove this 
conjecture, one should check that increasing the lattice size one gets smaller and smaller 
values of the mass gap. To get a clean indication, for the same value k. = 0.076, this will 
require the non-trivial task to compute the connected two-point correlator, with the same 
Swendsen-Wang cluster algorithm and a statistics of ~ 6 Msweeps, on a substantially 
larger lattice, say on a 48^. 

Before ending this section, 1 observe that from the results of Sect. 2 the non-perturbative 
infrared sector of the broken phase emerges as a threshold phenomenon starting at a mo- 
mentum scale 5 <^ M^. As anticipated, a massive free-field limit for the continuum theory 
implies a hierarchical structure of scales where the ratio ^ — ^ ^ ^ when t — ^ ^ oo. 
Therefore, the discrepancy in the values of the mass gap, that is already observed at 
K, = 0.076 by simply increasing L from 20 to 32, would show up on larger and larger 
lattice sizes L > 1/5 for k approaching the critical value Kc ~ 0.07484 [29]. 

For instance assuming a hierarchical relation of the type in Eq.(l), the minimum lattice 
size would be L > A/M^. Thus, using the relation [36] A ~ to express the ultraviolet 
cutoff of a A$4 theory in units of the inverse lattice spacing, one should expect 

L I.Stt , , 

-a > (W 

This gives f > 29, 52, 118, 471 for M/,a=0.4, 0.3, 0.2 , 0.1 respectively. In this sense, 
the range close to k, — 0.076, where M/jO ~0.4, is uniquely singled out since it hes in the 
scaling region but, at the same time, the required lattice sizes are not too large. For the 
same reason, the stability of a mass gap such as 0.205 when increasing the lattice size 
from = 36 to - = 48 (see Table 4 of Ref.[30]) might be consistent with the results 
presented here. In this case, in fact, to detect the same discrepancy found for M^a ~ 0.4 
a minimum value - 112 might be needed. 



4. Summary and conclusions 

In this paper, following the original idea of Ref.[3], I have presented additional evidences 
that the conventional perturbative picture of SSB might miss an important infrared phe- 
nomenon: the standard singlet Higgs boson might not be a purely massive field. In fact, 
as reviewed in the Introduction, there are arguments to expect a non-perturbative infrared 
sector corresponding to the long-wavelength excitations of the scalar condensate. 

On a more formal ground, as a consequence of the convexification process of the 
effective potential, the inverse zero-momentum connected propagator should be considered 
a two-valued function that, besides the standard massive solution G~^{p = 0) = M|, 
includes the value G~^^{p = 0) = as in a gapless theory. As shown in Ref. [12] and 
illustrated in Sect. 2, the convexification process is a threshold phenomenon that starts 
when the infrared cutoff k is below a tiny scale 5. In this regime, where the effective 
potential T4(0) deviates from the smooth semiclassical form of perturbation theory, the 
field strength starts to exhibit large deviations from unity with a strong peaking 

phenomenon that extends to the boundary of the fiatness region in the limit A; ^ 0. Such 
a behaviour could hardly be explained if the infrared region were that of a simple massive 
theory. 

In Sect. 3 (and in the Tables at the end of the paper) I have produced detailed 
numerical results from a lattice simulation of the broken phase in the 4D Ising limit of 
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the theory. They point to the following conclusions: 

i) differently from the symmetric phase, the single-particle energy spectrum is not 
well reproduced by (the lattice version of) the standard massive form -i^p^ + const, in the 
limit p — > and a |p|-independent mass parameter is only found at higher momenta 

ii) for the value k — 0.076, the mass gap reported by Balog et al.[30] on a 20^ lattice 
is mTs(O) = 0.392(1). On the other hand increasing the lattice size up to 32^ the data 
provide mTs(O) = 0.366(5), contrary to the expectation that there should be no significant 
change. At the same time, since the susceptibility is practically unchanged, this might 
represent one more evidence for the existence of two solutions for the zero-momentum 
propagator and for the subtle nature of the zero-momentum limit in the broken-symmetry 
phase. 

For this reason, the observed volume dependence of mTs(O), with its possible interpre- 
tation in terms of ('non-Goldstone') collective excitations of the scalar condensate, poses 
interesting questions and would deserve to be checked by other groups with new system- 
atic investigations of the L ^ oo hmit of the lattice theory. To this end, one should 
also take into account that the deviations from a simple massive spectrum a/p^ + 
might be confined to a region of momenta |p| < 6 where 6 vanishes in units of in 
the infinite-cutoff limit. This means that, approaching the continuum limit of the lattice 
theory, the volume dependence of the mass gap might require larger and larger lattice 
sizes before to show up. 

Acknowledgements 

I wish to thank Paolo Cea and Leonardo Cosmai for giving me all necessary ingredients 
for the lattice simulations. I also thank Dario Zappala for many useful discussions. 



12 



References 

[1] G. 't Hooft, In Search of the Ultimate Building Blocks, Cambridge University Press, 
1997, p.70. 

[2] M. ConsoU and P.M. Stevenson, Int. J. Mod. Phys. A15 (2000) 133, hep-ph/9905427. 

[3] M. Consoh, Phys. Lett. B512 (2001) 335; Phys. Rev. D65 (2002) 105217. 

[4] P. M. Stevenson, in Proceedings of the Conference on Violations of CPT and Lorentz 
invariance, Indiana 2001 (World Scientific, Singapore 2002). 

[5] G. E. Volovik, Phys. Rep. 351 (2001) 195. 

[6] M. Consoh, Phys. Lett. B541 (2002) 307. 

[7] G. E. Volovik, JETP Lett.73 (2001) 162. 

[8] F. Ferrer and J. Grifols, Phys. Rev. D63 (2001) 025020. 

[9] K. Symanzik, Commun. Math. Phys. 16 (1970) 48. 

[10] M. Consoh and E. Costanzo, Eur. Phys. J. C33 (2004) 297. 

[11] R. Jackiw, S. -Y. Pi and A. Polychronakos, Ann. Phys. 301 (2002) 157. 

[12] M. Consoli and D. Zappala, arXiv:hep-th/0606010, to appear in Phys. Lett. B. 

[13] F. J. Wegner and A. Houghton, Phys. Rev A8 (1973) 401. 

[14] T.S. Chang, D. D. Vvedensky and J.F. NicoU, Phys. Rep. 217 (1992) 280. 

[15] J. Polchinski, Nucl. Phys. B 231 (1984) 269. 

[16] C. Wetterich, Nucl. Phys. B352, 529 (1991); Phys Lett. B301 (1993) 90. 

[17] T. Morris, Int. J. Mod. Phys. A 9 (1994) 2411; Phys. Lett. B329 (1994) 241. 

[18] J. Berges, N. Tetradis and C. Wetterich, Phys.Rep 363 (2002) 223. 

[19] K.-I. Aoki, A. Horikoshi, M. Taniguchi and H. Terao, The Exact Renormalization 
Group, Proceedings of the Workshop on the Exact Renormalization Group, Faro, 
Portugal, September 10-12 1998, (World Scientific 1999), arXiv:hep-th/9812050; A. 
S. Kapoyannis and N. Tetradis, Phys. Lett. A 276 (2000) 225. 

[20] D. F. Litim, J. M. Pawlowski and L. Vergara, Convexity of the effective action from 
functional flows , arXiv:hep-th/0602140. 

[21] J. Alexandre, V. Branchina and J. Polonyi, Phys. Lett. 445 (1999) 351. 

[22] D. Zappala, Phys. Lett. A290 (2001) 35. 

[23] A. Bonanno and D. Zappala, Phys. Lett. B504 (2001) 181; M.Mazza and D. Zappala, 
Phys. Rev. D64 (2001) 105013. 



13 



[24] D.F. Litim and J. M. Pawlowski, Phys. Lett. B546 (2002) 279; Phys. Rev. D66 
(2002) 025030. 

[25] D.F. Litim and J. M. Pawlowski, Phys. Lett. B516 (2001) 197; Phys. Rev. D65 
(2002) 081701. 

[26] For a review of the rigorous results, see R. Fernandez, J. Frohhch, and A. D. 
Sokal, Random Walks, Critical Phenomena and Triviality in Quantum Field The- 
ory ( Springer- Verlag Berlin Heidelberg 1992). 

[27] P. Cea. M. Consoh, L. Cosmai and P. M. Stevenson, Mod. Phys. Lett. A14 (1999) 
1673, hep-lat/9902020. 

[28] P. M. Stevenson, Nucl. Phys. B729 (2005) 542. 

[29] D.S. Gaunt, M.F. Sykes and S. McKenzie, J. Phys. A12 (1979) 871; D. Stauffer and 
J. Adler, Int. J. Mod. Phys. C8 (1997) 263. 

[30] J. Balog et al, Nucl. Phys. B714 (2005) 256, hep-lat/0412015. 

[31] I. Montvay and G. Miinster, Quantum Fields on a Lattice, Cambridge University 
Press, 1994. 

[32] R.H. Swendsen and J.S. Wang, Phys. Rev. Lett. 58 (1987) 86. 

[33] G. Frick, K. Jansen, and P. Seuferhng, Phys. Rev. Lett. 63 (1989) 2613. 

[34] B. Efron, Jackknife, the Bootstrap and other resampling plans (SIAM, 1982); B. A. 
Beg and A. H. Billoire, Phys. Rev. D 40 (1989) 550. 

[35] J. Balog, F. Niedermayer and P. Weisz, hep-lat/0601016. 

[36] D. E. Brahm, arXiv:hep-lat/9403021. 



14 



t Ci{t, 0; k = 0) statistical error 



0.46445456707497D-03 

1 0.37571910789012D-03 

2 0.30382577283847D-03 

3 0.24582066021523D-03 

4 0.19905074016424D-03 

5 0.16136544145195D-03 

6 0.13104967814461D-03 

7 0.10671304049363D-03 

8 0.87238915133100D-04 

9 0.71738763750074D-04 

10 0.59502097899827D-04 

11 0.49970839560653D-04 

12 0.42700968213871D-04 

13 0.37366513141293D-04 

14 0.33719361396463D-04 

15 0.31592799783943D-04 

16 0.30893472344528D-04 

17 0.31592799783943D-04 

18 0.33719361396463D-04 

19 0.37366513141293D-04 

20 0.42700968213871D-04 

21 0.49970839560653D-04 

22 0.59502097899827D-04 

23 0.71738763750074D-04 

24 0.87238915133100D-04 

25 0.10671304049363D-03 

26 0.13104967814461D-03 

27 0.16136544145195D-03 

28 0.19905074016424D-03 

29 0.24582066021523D-03 

30 0.30382577283847D-03 

31 0.37571910789012D-03 

32 0.46445456707497D-03 



0.46037664989336D-06 

0.46116594635711D-06 
0.45018527889238D-06 
0.43438332218099D-06 
0.41560747498003D-06 
0.39542065913881D-06 
0.37531332727689D-06 
0.35557753750683D-06 
0.33706071111197D-06 
0.32019883688822D-06 
0.30574726422856D-06 
0.29433390224131D-06 
0.29772329849853D-06 
0.29930693571779D-06 
0.29930065059117D-06 
0.29950600189874D-06 
0.30005725246790D-06 
0.29950600189874D-06 
0.29930065059117D-06 
0.29930693571779D-06 
0.29772329849853D-06 
0.29433390224131D-06 
0.30574726422856D-06 
0.32019883688822D-06 
0.33706071111197D-06 
0.35557753750683D-06 
0.37531332727689D-06 
0.39542065913881D-06 
0.41560747498003D-06 
0.43438332218099D-06 
0.45018527889238D-06 
0.46116594635711D-06 
0.46037664989336D-06 



t — range x^/d.o.f. £'(k) mTs(k) 

0^32 0.072 0.000 0.2129 (2) 0.2133 (2) 

1- 31 0.075 0.000 0.2128 (2) 0.2133 (2) 

2- 30 0.052 0.000 0.2128 (3) 0.2132 (3) 

3- 29 0.035 0.000 0.2127 (3) 0.2131 (3) 

Table 1: The raw correlator data obtained for a spatial momentum k = 0, on a 32^ lattice, 
for the value k, = 0.074 in the symmetric phase. The statistics is 135Ksweeps. I also show 
the quality of the fits, the effective single-particle energy and the time-slice mass obtained 
by using Eqs.(29) and (26) in various ranges of t. 
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_t Ci(t.O:k) 

0.14623108565057D-03 

1 0.77214443895274D-04 

2 0.40796525737581D-04 

3 0.21558338760062D-04 

4 0.11389647365225D-04 

5 0.60169117411735D-05 

6 0.31785227732011D-05 

7 0.16760427054163D-05 

8 0.88228124131541D-06 

9 0.46368972631045D-06 

10 0.24302326620048D-06 

11 0.12847759060431D-06 

12 0.68535903339602D-07 

13 0.37768265153790D-07 

14 0.22751593190996D-07 

15 0.15386378604173D-07 

16 0.13380756871513D-07 

17 0.15386378604173D-07 

18 0.22751593190996D-07 

19 0.37768265153790D-07 

20 0.68535903339602D-07 

21 0.12847759060431D-06 

22 0.24302326620048D-06 

23 0.46368972631045D-06 

24 0.88228124131541D-06 

25 0.16760427054163D-05 

26 0.31785227732011D-05 

27 0.60169117411735D-05 

28 0.11389647365225D-04 

29 0.21558338760062D-04 

30 0.40796525737581D-04 

31 0.77214443895274D-04 

32 0.14623108565057D-03 



statistical error 
0.14165730212821D-07 
0.1328835521 1454D-07 
0.11455026260932D-07 
0.95795734647824D-08 
0.81678791276023D-08 
0.78549959762898D-08 
0.64601886340349D-08 
0.57641678003130D-08 
0.53085501360843D-08 
0.46940187418747D-08 
0.42367563032742D-08 
0.38866243669447D-08 
0.37374075184433D-08 
0.37369789004914D-08 
0.39349497056091D-08 
0.42589381877184D-08 
0.44421054740855D-08 
0.42589381877184D-08 
0.39349497056091D-08 
0.37369789004914D-08 
0.37374075184433D-08 
0.38866243669447D-08 
0.42367563032742D-08 
0.46940187418747D-08 
0.53085501360843D-08 
0.57641678003130D-08 
0.64601886340349D-08 
0.78549959762898D-08 
0.81678791276023D-08 
0.95795734647824D-08 
0.11455026260932D-07 
0.13288355211454D-07 
0.14165730212821D-07 



t — range x^/d.o.f. mTs(k) 

(r32 0.589 0.375 0.6382 (1) 0.2142 (3) 

1- 31 0.377 0.375 0.6380 (2) 0.2135 (5) 

2- 30 0.382 0.375 0.6381 (3) 0.2139 (8) 

3- 29 0.328 0.375 0.6384 (4) 0.2149 (13) 

Table 2: The raw correlator data obtained for a spatial momentum k^ = 0.375 (corre- 
sponding to the integer assignment {nx,ny,nz) = (0,3,1) on a 32^ lattice), for the value 
K, — 0.074 in the symmetric phase. The statistics is 135Ksweeps. I also show the quahty 
of the fits, the effective single-particle energy and the time-slice mass obtained by using 
Eqs. (29) and (26) in various ranges of t. 
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t 




Ci(t,0;k) 




statistical error 







0. 


,37264797340762D-03 


0, 


.39868629339350D- 


-07 


1 


0. 


, 144237037531 17D-03 


0, 


.35128421137473D- 


-07 


2 


0. 


,55852686582681D-04 


0, 


.29204927503163D- 


■07 


3 


0. 


,21643050470909D-04 


0, 


.25324204561198D- 


-07 


4 


0, 


,83911901659690D-05 


0, 


.23695231347215D- 


-07 


5 


0. 


,32553421665971D-05 


0, 


.23752227353389D- 


-07 


6 


0. 


,12625860175981D-05 


0, 


.25644772176930D- 


-07 


7 


0. 


,48959050096638D-06 


0, 


.23531686001271D- 


■07 


8 


0. 


,19263769002945D-06 


0, 


.16045711212737D- 


-07 


9 


0. 


,69459096205410D-07 


0, 


.17742752632037D- 


-07 


10 


0. 


,30668878114950D-07 


0, 


.19019731595239D- 


-07 


11 


0. 


,69459099077044D-07 


0, 


.17742753057175D- 


-07 


12 


0. 


,19263769455122D-06 


0, 


.16045711766908D- 


■07 


13 


0, 


,48959049739535D-06 


0, 


.23531683326186D- 


-07 


14 


0. 


,12625860091612D-05 


0, 


.256447747460 16D- 


-07 


15 


0. 


,32553421642504D-05 





.23752229326030D- 


-07 


16 


0. 


,83911901631759D-05 


0, 


.23695234653985D- 


-07 


17 


0. 


,21643050344165D-04 


0, 


.25324219112844D- 


■07 


18 


0. 


,558526864859360-04 


0, 


.29204909466865D- 


-07 


19 


0, 


,14423703585060D-03 





.35128469121253D- 


-07 


20 


0. 


,37264797340762D-03 


0, 


.39868629339350D- 


-07 



— range 


xVd-O.f. 






mTs(k) 


0-20 


0.451 


0.922 


0.9489 (1) 


0.2186 (6) 


1-19 


0.274 


0.922 


0.9485 (3) 


0.2162 (14) 


2-18 


0.244 


0.922 


0.9479 (6) 


0.2133 (33) 


3-17 


0.282 


0.922 


0.9476 (16) 


0.2120 (83) 



Table 3: The raw correlator data obtained for a spatial momentum = 0.922 (corre- 
sponding to the integer assignment {nx,ny,nz) = (3,1,0) on a 20^ lattice), for the value 
K, — 0.074 in the symmetric phase. The statistics is 215 Ksweeps. I also show the quahty 
of the fits, the effective single-particle energy and the time-slice mass obtained by using 
Eqs.(29) and (26) in various ranges of t. 
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t Ci{t, 0; k = 0) statistical error 



0.87514834087643D-03 

1 0.61452112729030D-03 

2 0.41664497529767D-03 

3 0.28249401086311D-03 

4 0.19223970120763D-03 

5 0.13173539366874D-03 

6 0.91507623807116D-04 

7 0.65276805609368D-04 

8 0.48943304332497D-04 

9 0.40073467468121D-04 

10 0.37300092795145D-04 

11 0.40073467468121D-04 

12 0.48943304332497D-04 

13 0.65276805609368D-04 

14 0.91507623807116D-04 

15 0.13173539366874D-03 

16 0.19223970120763D-03 

17 0.28249401086311D-03 

18 0.41664497529767D-03 

19 0.61452112729030D-03 

20 0.87514834087643D-03 



0.40185555074734D-05 
0.39930724185063D-05 
0.39573238547867D-05 
0.39576325347383D-05 
0.39509605466453D-05 
0.39095819880648D-05 
0.39253975622550D-05 
0.39656063145965D-05 
0.40437233918344D-05 
0.40949224665482D-05 
0.41228418720897D-05 
0.40949224665482D-05 
0.40437233918344D-05 
0.39656063145965D-05 
0.39253975622550D-05 
0.39095819880648D-05 
0.39509605466453D-05 
0.39576325347383D-05 
0.39573238547867D-05 
0.39930724185063D-05 
0.40185555074734D-05 



t — range x^/d.o.f. E(k) mTs(k) 

0^20 r572 0.000 0.3804 (17) 0.3827 (17) 

1- 19 0.005 0.000 0.3895 (24) 0.3920 (24) 

2- 18 0.005 0.000 0.3892 (35) 0.3917 (35) 

3- 17 0.002 0.000 0.3884 (50) 0.3909 (51) 

Table 4: The raw correlator data, at k = and for k, — 0.076 in the broken phase. The 
statistics is 6 Msweeps on a 20^ lattice. I also show the quality of the fits with Eq.(29), the 
effective single particle energy and the time-slice mass obtained in various ranges of t. The 
vacuum expectation value was = 0.30157(3) and the susceptibility xiatt = 37.87(9) . 
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E{k) 


mTs(k) 
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0.4024 (13) 
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0.098 


0.5019 (16) 


0.3992 (21) 
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0.098 


0.5014 (25) 


0.3985 (32) 
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0.5014 (38) 


0.3985 (49) 



Table 5: The raw correlator data for k, — 0.076 in the broken phase, on a 20^ lattice, for 
the integer assignment (1,0,0). The statistics is 230 Ksweeps. I also show the quality of 
the fits, the effective single-particle energy and the time-slice mass obtained in various 
ranges of t. 
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Ci(t,0;k) 



statistical error 



0.60134431492201D-03 

1 0.33121256195148D-03 

2 0.18340753552086D-03 

3 0.10169038556086D-03 

4 0.56548852939487D-04 

5 0.31288412716245D-04 

6 0.17273284231833D-04 

7 0.98445753758052D-05 

8 0.58757929618968D-05 

9 0.40294379323518D-05 

10 0.34202276827878D-05 

11 0.40294380412313D-05 

12 0.58757929936507D-05 

13 0.98445753579537D-05 

14 0.17273284284194D-04 

15 0.31288412838347D-04 

16 0.56548852916696D-04 

17 0.10169038565092D-03 

18 0.18340753431478D-03 

19 0.33121256136276D-03 

20 0.60134431492201D-03 



0.72278071460665D-06 
0.67224614280968D-06 
0.57001952486459D-06 
0.58126295776147D-06 
0.62601000175195D-06 
0.61885213006675D-06 
0.61847800661692D-06 
0.60619584062421D-06 
0.64071878068011D-06 
0.68994571329391D-06 
0.74746648254179D-06 
0.68994574565471D-06 
0.64071876223713D-06 
0.60619580358585D-06 
0.61847800602069D-06 
0.61885209771265D-06 
0.62600984641147D-06 
0.58126288128429D-06 
0.57001966754573D-06 
0.67224582018553D-06 
0.72278071460665D-06 



t — range x^/d.o.f. -^(k) mTs(k) 

0^20 (1358 0196 0.5932 (8) 0.4081 (12) 

1- 19 0.036 0.196 0.5905 (14) 0.4039 (21) 

2- 18 0.036 0.196 0.5899 (23) 0.4030 (35) 

3- 17 0.042 0.196 0.5902 (42) 0.4034 (63) 

Table 6: The raw correlator data for k, — 0.076 in the broken phase, on a 20^ lattice, for 
the integer assignment (1,1,0). The statistics is 230 Ksweeps. I also show the quality of 
the fits, the effective single-particle energy and the time-slice mass obtained in various 
ranges of t. 
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0.4163 (6) 
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0.278 
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0.4122 (15) 
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0.272 
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0.9609 (15) 


0.4150 (40) 



Table 7: The raw correlator data for k, — 0.076 in the broken phase, on a 20^ lattice, for 
the integer assignment (3,0,0). The statistics is 10 Msweeps. I also show the quality of 
the fits, the effective single-particle energy and the time-slice mass obtained in various 
ranges of t. 
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Table 8: The raw correlator data for k, — 0.076 in the broken phase, on a 20^ lattice, for 
the integer assignment (3,1,0). The statistics is 10 Msweeps. I also show the quality of 
the fits, the effective single-particle energy and the time-slice mass obtained in various 
ranges of t. 
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Table 9: The raw correlator data for k — 0.076 in the broken phase, on a 20^ lattice, for 
the integer assignment (3,1,1). The statistics is 10 Msweeps. I also show the quality of 
the fits, the effective single-particle energy and the time-slice mass obtained in various 
ranges of t. 
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Table 10: The raw correlator data at k = obtained from a simulation on a 32"^ lattice, for 
K = 0.076 in the broken phase. The statistics is 3500 Ksweeps. Only data with S/N > 1 
are reported. The vacuum expectation value was = 0.30158(2) and the susceptibility 
Xiatt = 37.70(11). 
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Table 11: The raw correlator data at k = obtained from a simulation on a 32"^ lattice, for 
K = 0.076 in the broken phase. The statistics is 1250Ksweeps. Only data with S/N > 1 
are reported. The vacuum expectation value was = 0.30156(4) and the susceptibility 
Xiatt = 37.72(18). 
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Table 12: The raw correlator data at k = obtained from a simulation on a 32'^ lattice, 
for K, — 0.076 in the broken phase. The statistics is 820Ksweeps. Only data with S/N > 1 
are reported. The vacuum expectation value was = 0.30157(3) and the susceptibihty 
Xiatt = 37.74(20). 
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Table 13: The raw correlator data at k = from a simulation on a 32^ lattice, for 
K — 0.076 in the broken phase. The statistics is 620Kswccps. Only data with S/N > 1 
are reported. The vacuum expectation value was — 0.30157(3) and the susceptibility 
Xiatt = 37.72(28). 



t — range x^/d-o.f. k^ -^(k) mTs(k) 

0^32 0.779 0.000 0.3443 (25) 0.3460 (25) 

1- 31 0.255 0.000 0.3624 (37) 0.3644 (37) 

2- 30 0.254 0.000 0.3655 (54) 0.3675 (54) 

3- 29 0.256 0.000 0.3622 (77) 0.3642 (77) 

4- 28 0.252 0.000 0.3546 (110) 0.3565 (110) 

Table 14: The results of the global fit to the correlator data reported in Tables 10-13, 
for K = 0.076 in the broken phase. The total statistics is 6.19 Msweeps on 32^ lattices. 
I report the quality of the fits with Eq.(29), the effective single particle energy and the 
time-shce mass obtained in various ranges of t. These results should be compared with 
those in Table 4 from the 20^ lattice. 
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t — range x^/d-o.f. -E'i(k) mTs(k) 

0^20 1.573 0.000 0.3804 (17) 0.3827 (17) 

1- 19 0.002 0.000 0.3880 (71) 0.3904 (71) 

2- 18 0.001 0.000 0.3866 (113) 0.3890 (113) 

3- 17 0.000 0.000 0.3861 (159) 0.3885 (159) 

Table 15: The results of a fit to the correlator data reported in Table 4 from the 20^ 
lattice. The fit is performed using a constrained 2- mass formula with £'2(0) = 2£'i(0) 
as suggested in Ref.[35]. I report the quality of the fits, £"1(0) and the time-slice mass 
obtained in various ranges of t. These results should be compared with those in Table 4 
obtained fitting the correlator data to Eq.(29). 



t — range x^/d.o.f. k^ -^i(k) ^Ts(k) 

0^32 0?787 000 0.3443 (25) 0.3460 (25) 

1- 31 0.258 0.000 0.3624 (37) 0.3644 (37) 

2- 30 0.236 0.000 0.3416 (198) 0.3433 (198) 

3- 29 0.225 0.000 0.3143 (324) 0.3156 (324) 

Table 16: The results of a fit to the correlator data reported in Tables 10-13 from the 
32^ lattice. The fit is performed using a constrained 2-mass formula with £2(0) = 2£'i(0) 
as suggested in Ref.[35]. 1 report the quality of the fits, £1(0) and the time-slice mass 
obtained in various ranges of t. These results should be compared with those in Table 
14 obtained fitting the same correlator data with Eq.(29) and with those in Table 15 
obtained from the 20^ lattice data using the same 2-mass fitting function. 
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